Abstract 



In this paper we present the results of computer searches using 
a variation of an energy minimization algorithm used by Kottwitz 
for finding good spherical codes. We prove that exact codes exist 
by representing the inner products between the vectors as algebraic 
numbers. For selected interesting cases, we include detailed discussion 
of the configurations. Of particular interest are the 20-point code in 
M 6 and the 24-point code in M. 7 , which are both the union of two 
cross polytopes in parallel hyperplanes. Finally, we catalogue all of 
the codes we have found. 
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1 Introduction 

Given N points that lie on the unit sphere S 1 ™ -1 in M n , we wish to determine 
how they should be placed so that the minimal distance between any two 
points is maximized. Any set of points on the unit sphere is called a spherical 
code, and the problem of finding the best code has been proposed many 
times and in many contexts (though usually only in three dimensions), from 
packing circles on a sphere to distributing orifices on pollen-grains (e.g. [5], 
|18j). The only known optimal solutions in three dimensions are for N < 12 
and iV = 24, while in four dimensions only the cases iV < 8, N = 10 pQ, 
and iV = 120 have been proven. The remaining known optimal codes are for 
N < 2n for any n, in which the points form either a simplex or a subset of 
the cross polytope; certain codes derived from the Leech lattice and the E$ 
root system; and lastly an infinite family based on isotropic subspaces [2]. 
Excluding these few cases, the best known codes have been found either by 
specific constructions or, more commonly, by computer searches using various 
optimization algorithms. Given the difficulty of proving the optimality of 
even very small codes, most work related to this problem has been in finding 
close approximations to good configurations. The most extensive table of 
codes was created by N. J. A. Sloane, with the collaboration of R. H. Hardin, 
W. D. Smith and others, and is available electronically |17j . 

In this work we use the technique of energy minimization to find good 
spherical codes. Leech first observed the possibility of such an approach, 
and it has been used in previous work several times [T5] . Kottwitz gave a 
fairly comprehensive list of three-dimensional codes for up to 90 points [T2] , 
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which was expanded later by Buddenhagen and Kottwitz while searching for 
codes with multiplicity greater than one (i.e. codes for which there exist 
distinct configurations that obtain the same optimal distance) jl]. Nurmela 
investigated some global optimization methods based on energy minimiza- 
tion and provided numerical results for codes in up to five dimensions [15j, 
but the most successful implementations have been based on using a large 
number of random starts and a local optimization algorithm so that there 
is a high probability that at least one of them converges to the global opti- 
mum. Our algorithm has several changes to improve on old implementations 
as well, such as reducing exponent bias and choosing a different local opti- 
mization method. This has resulted in three improved spherical codes in four 
dimensions as well as new higher dimensional codes, particularly in six and 
seven dimensions, that exhibit interesting symmetries. Section [2] describes 
our algorithm in detail and Section [3] gives a brief analysis of the improved 
codes. 

The remainder of this work focuses on providing rigorous analysis of 
codes, particularly showing that the best known codes can be represented 
as exact configurations in terms of algebraic numbers as opposed to close 
estimates. This is an important step toward rigorizing computer solutions to 
the problem and allows us to observe true equalities and relations between 
points and edges in the code. Buddenhagen and Kottwitz did similar work 
while looking for three-dimensional codes with multiplicity greater than one, 
and provided detailed exact descriptions of the two distinct optimal 15-point 
codes jl]. In several cases we were able to identify a considerable amount 
of underlying symmetry and structure using the methods of Section HI For 
three particularly interesting cases we provide brief discussion of the config- 
urations in Section The final section gives tables of the exact codes, in 
dimensions four through eight, based on their minimal polynomials. 



We consider an inverse power law potential function on the spherical code 
C = {xi, x 2 , . . . , a; at} C S"" -1 defined by 



2 Methodology 




(i) 
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where a is a constant to prevent overflow, s is the exponent of the inverse 
power law, and \x{ — Xj\ is the Euclidean distance between the two points 
in W 1 . As s — > oo, the smallest distance will be the dominating term in the 
energy expression. Hence, minimizing this potential function over all possible 
C will give approximate codes that get better as s increases. In certain highly 
symmetrical cases, when s is sufficiently large, the method will produce the 
actual optimal code, whereas in general it converges to the optimal code as 
s — > oo. 

Several issues arise, however, if one tries to minimize E using a large s 
with a random configuration of points. There is an implementation problem, 
in that a cannot be chosen well so that E neither overflows nor is too small. 
It has no abstract mathematical role, but it is important to consider when 
dealing with floating-point arithmetic. If a poor choice of a is used with a 
power of s — 1, 000, 000 for example, a slight deviation in \xj —Xj \ from a by 
only 0.001 will result in either 1.001 1000000 « 10 434 or 0.999 1000000 « 10" 435 , 
neither of which would be handled well by a computer. There is also a 
problem regarding the behavior of the function E as s increases because the 
number of local minima increases as well. Starting with a high exponent in 
the inverse power interaction will usually converge quickly to a poor local 
minimum. 

Taking into consideration these problems, Kottwitz used the following 
algorithm [T2] : 

1. Start with s = 80 and a random configuration of points; 

2. Run a local optimization algorithm until convergence; 

3. Double s and repeat the previous step starting from the local 
optimum found, stopping after the optimization is run for s = 
1,310,720. 

This remedies the issues described because a can be chosen before each opti- 
mization to reflect the minimal distance of the configuration, and the points 
will tend toward a good code without getting stuck early on in a bad lo- 
cal minimum. For our work we use a similar strategy, but to avoid the bias 
introduced by fixing an initial exponent we start with a random exponent be- 
tween 10 and 160 each time. Furthermore, we allow the doubling procedure 
to continue until the exponent exceeds 100, 000, 000. Before each optimiza- 
tion procedure, we choose a to be the exact minimal distance between all 
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pairs of points, as done by Nurmela [JS]. This is a good choice because as 
the power increases the improvement in the minimal distance decreases, and 
we have found that this balance keeps the energy expression in the desired 
range. 

While Kottwitz used gradient descent and Newton-Raphson methods to 
optimize the energy function at each step jl], [12], we picked the nonlin- 
ear conjugate gradient algorithm for several reasons. Firstly, it is easy to 
implement yet much more effective than gradient descent. It also only re- 
quires the gradient at each step and does not need to compute the Hessian, 
which is costly and very complex especially as the dimension increases. We 
experimented with three versions of the nonlinear conjugate gradient algo- 
rithm: Fletcher-Reeves, Polak-Ribiere, and Hestenes-Stiefel. We found that 
the Hestenes-Stiefel formula for updating the conjugate direction was most 
effective for our model, converging faster than Polak-Ribiere and displaying 
better results than Fletcher-Reeves. 

We must also consider the constrained nature of optimization on a sphere, 
as the conjugate gradient algorithm does not accommodate for constraints on 
the domain. Our solution to this is to, after each movement in the conjugate 
direction, simply scale each point in C back onto the sphere. We realize, 
however, that as the code approaches a local minimum, the points will tend 
to move almost directly away from the sphere because most of the movement 
along the tangent space will be balanced by their neighbors. When scaled 
they will consequently appear to have moved very little if at all, slowing con- 
vergence and possibly even getting stuck before converging. To compensate 
for this, when the exponent in the inverse power law interaction is high, for 
each point we only look at the component of the gradient that lies in the 
tangent space to the sphere at that point. We do not implement such a 
procedure for lower exponents because the conjugate directions depend on 
previous gradients as well so the magnitudes should be consistent relative to 
each other. 

Using the method that has been described, implemented in C++, we ran 
the optimization with at least 1000 random starts. Taking the best result 
of these runs, we then applied a similar optimization to the configuration 
to improve the precision using a starting exponent of 10, 000, 000 and dou- 
bling until 640, 000, 000. Once this converged we were able to obtain fairly 
accurate results for our codes, but not accurate enough to find their exact 
representations. Finally, we found the set of shortest edges to approximately 
five or six decimal places and set them equal to each other, which determined 
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a rigid structure solved using Newton-Raphson methods for a small number 
of iterations to a precision of about 500 decimal digits. Then we placed the 
remaining points, the rattlers, by moving them iteratively as far as possible 
from the rigid points, again maximizing the minimal distance. We used the 
PARI/GP software for the last two optimization steps. 

After the full procedure, our set of codes performed very well against 
previous data. For the three-dimensional codes we ran our program on, we 
reproduced the best known results in every case, while for four-dimensional 
codes we were able to find three codes with 40, 68, and 71 points that im- 
proved on the best previous results. In higher dimensions for smaller numbers 
of points we are confident that we were able to find good (and most likely 
optimal) configurations. 

3 Improved Codes in M 4 

Our work has improved on three codes in four dimensions, namely those with 
40, 68, and 71 points. They are briefly examined in this section. Note that 
the improvements are small - on the order of l/1000th or less. A catalogue 
of all of the codes can be found online at 

http : //www . aimath . org/ data/paper/WangSphericalCode/. 

The data available online reflects only results prior to the optimization done 
in PARI/GP. 

3.1 40 points 

The maximal cosine of 0.65049780106271773133. . . occurs 108 times, but is 
the only inner product that occurs more than twice, suggesting that this 
configuration is quite asymmetric. 

3.2 68 points 

The maximal cosine of 0.75104449257228207352... occurs 196 times, and 
there is one other inner product - 0.74925795575260304153 ... - that occurs 
three times. Every other edge length appears to be unique. 
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3.3 71 points 



The maximal cosine of 0.75637601134814761871 . . . occurs 199 times, and 
there are three other inner products that occur three times each in the 
code; they are 0.75152318440020477595 . . ., 0.75574916415250115097 . . ., and 
0.75632115855377282465 . . .. The other edges are unique. 

4 Properties for Identifying Codes 

We examine some properties of spherical codes that can hopefully character- 
ize them more substantially and rigorously than a list of floating-point coordi- 
nates. First, we introduce the NxN Gram matrix G such that G^ = Xj), 
where (x, y) denotes the standard inner product between the unit vectors x 
and y. Since all vectors on the unit sphere have magnitude one, we also have 
Gij = cos (9ij), where 6ij is the angle between the vectors Xi and Xj. 

4.1 Edge Length Occurrences 

One effective way of determining whether or not a spherical code has sym- 
metry is to look at the number of occurrences of each edge length or the 
number of unique edge lengths between distinct points, which we will call m. 
Since the distance and inner product are related by 

\xi Xj | =2 2 (iCj, Xj) , (2) 

this consists of simply looking at the unique entries of the Gram matrix 
and the number of times each is repeated. We note that considering edge 
length occurrences is a generic and flexible strategy for describing codes and 
detecting overall symmetry, especially when m is small relatively to the code 
size. To find more specific symmetries, we can look at the edge lengths 
emanating from any one point in the code and check if any other points 
shared these. In particular, we can often group vertices together as equivalent 
if the lists of edge lengths emanating from them are exactly the same. We 
note that this is a necessary but not sufficient condition for the vertices to 
be equivalent under the action of the symmetry group of the configuration. 
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4.2 Exact Algebraic Numbers 

One inevitable concern that arises from computer searches for spherical codes 
is whether or not the set of points found by the computer really exists. 
Numerical results can be useful for recognizing the exact solutions, but with 
finite precision it is not clear if all the distances are really what we think they 
are. Consequently it is important to have some sort of method of converting 
the work of computer approximations into definitive arrangements so that 
they may be analyzed with more certainty and have hope of being proven 
optimal as well. 

Since a code is the solution to a number of polynomial equalities between 
the shortest edges, the coordinates of each rigid point in the code are algebraic 
and lie in a number field. Moreover, we then know that the inner products 
are algebraic as well. Once we have determined the coordinates of the points 
in the code to high precision we are able to effectively compute, using the 
PARI/GP math software command algdep, the minimal polynomial for the 
maximal inner product between two distinct vectors, i.e., 

u = max (xi, Xj) 

l<i<j<N X ' 31 

This was done for a number of three-dimensional codes by Buddenhagen and 
Kottwitz, and we extend their work. Just looking for a minimal polynomial 
P so that P{u) pa can be misleading, however, because even with high 
precision it could be the case that u is in fact not the root of P, but just 
very close. On the other hand, if we could show that the code with exactly 
the root of P (closest to u, to be specific) as the maximal inner product does 
exist, we would have an exact arrangement that is practically equivalent. We 
have been able to do so for most codes that we could calculate the minimal 
polynomial of. 

To verify the existence of the codes, we first use the PARI/GP command 
lindep to specify the elements of G in terms of polynomials in some primitive 
element of the number field. Empirically we have found that u is often 
sufficient for doing this, though in several cases this has not been true (in 
particular for 22 points in three dimensions, 21 points in five dimensions, and 
14 points in six dimensions). Once we have the exact Gram matrix, we may 
check if G is indeed a Gram matrix of a set of points on the unit sphere in 
W n , namely that it satisfies the following properties: (a) G is symmetric; (b) 
the diagonal elements of G are all equal to 1; (c) G has rank at most n; and 
(d) G is positive semidefinite. If G satisfies each of these conditions, it must 
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(3) 



be the Gram matrix of a spherical code, and thus we are able to prove the 
existence of such a code. The first two properties are easy to check, and the 
last two can be tested by observing the characteristic polynomial Pg(^) of 
G, which will be of degree N. Let 

Pg(A) = a N \ N + a N ^X N - 1 + • • • + a x X + a . (4) 

It is known that G has rank at most n if and only if pc(X) has a root at A = 
with multiplicity at least N — n, or equivalently a = a\ = ■ ■ ■ = a^^ n ^i = 0. 
Also, G is positive semidefinite if and only if the remaining coefficients are 
nonzero and alternate in sign. This makes it relatively easy to check if G is 
indeed a Gram matrix. 



5 Discussion of Selected Cases 

Here we present more detailed results of three cases which we found partic- 
ularly interesting and which have not been previously described in the lit- 
erature (except the last). Some brief discussion of other cases can be found 
online at http://www.aimath.org/data/paper/WangSphericalCode/. 

5.1 12 points in R 5 

This is a very simple configuration, so it is not surprising that the Gram 
matrix is quite simple as well, consisting of only five unique entries (excluding 
the l's along the diagonal). It is shown here: 
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where V and W satisfy V > W and are given in terms of u by 
V = \u 2 -u-\, W = -5u 2 -4u. 



There are two aspects of this Gram matrix that are particularly intrigu- 
ing. The first one that stands out is that all the edges between V\ and V 2 
correspond to an inner product of u or —u, which is a result of the points 
in V2 being antipodal. The second is that every element is represented as 
a polynomial in u of degree two or less. The minimal polynomial, however, 
is 25-u 4 + 30-u 3 + 24-u 2 + 2u — 1, which has degree four. This is unusual be- 
cause in most cases a significant number of elements of the Gram matrix 
are polynomials in u of maximal degree, i.e. one less than the degree of the 
minimal polynomial. While we do not have an explanation for this, it may 
be significant to the structure of the code. 

5.2 20 points in R 6 

This configuration, the first especially nice one, is highly symmetric; each 
point in the entire configuration is equivalent to every other one. Each point 
is at the minimal distance with 11 other points in the code. Furthermore, 
since the minimal polynomial is 14m — 3, all of the entries in the Gram 
matrix are rational numbers, specifically multiples of jj. The dot products 
corresponding to the edges that emanate from any point are (in order from 
largest to smallest) 

r 3 1 5349 

\ ' 14' '~14'"l4'"7'~7'~14 
3 

The only repetitions are — , which shows up 11 times, and — which shows 
up three times. 

It turns out that this code is the union of two five- dimensional cross 
polytopes in parallel hyperplanes. We can orient them with respect to one 
another: if we choose coordinates to allow one of the cross-polytopes to 
be given by plus or minus the standard orthonormal basis vectors, then 
after projection into the equatorial hyperplane and rescaling, the other cross 
polytope is given by 
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VTT' _ TT'~TT'TT'~TTy 

_ / 6 2 6 6 3 \ 

V IT' IT' IT' IT' — IT y 

/ 3 6 2 6 6 \ 

" V TT'~IT' IT'Ti'Tiy 

/ 6 3 6 2 6 \ 

\JT IT' - TT' IT' TTy 

/ 6 6 3 6 2 \ 

" \Ji' - Ti'Ti' - Ti'TiJ 



The way the orientation is determined is still unknown, but there is evi- 
dence of some sort of underlying structure that is worth investigating further. 



5.3 24 points in R 7 

This is the other particularly nice configuration. All 24 points in this code 
are symmetric with respect to one another and the minimal polynomial is 
19m 2 + 2u — 1, so the Gram matrix has very few entries and they are all 
simple. The only entries other than the l's along the diagonal are u, —u, 
—3u, and 2u — 1. Each vector has inner product u with 15 other vectors, 
— u with 2 other vectors, —3m with 5 other vectors, and 2u — 1 with a single 
other vector. This configuration appears to be a very good code as both the 
22- and 23-point codes are the same as this with points removed. 

It is very similar to the 20-point code in R e as it also consists of two cross- 
polytopes in parallel hyperplanes. As with the other case, we can orient the 
cross polytopes relative to each other by setting one to be plus or minus 
the standard orthonormal basis vectors. After projection into the equato- 
rial hyperplane and rescaling, the other cross polytope is then given by the 
following vectors: 
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This is a very structured set of vectors, but also differs from the vectors 
in the six dimensional case. It would be an interesting question to determine 
how in general the cross polytopes should be arranged relative to each other 
to maximize the minimal distance. 



6 Catalogue of Exact Configurations 

We present several tables of results for newly calculated exact spherical codes 
in dimensions four through eight, specifically the maximal inner product u 
and the minimal polynomial. In each dimension, we attempted to calculate 
a minimal polynomial for the maximum inner product of the configurations 
for up to 24 points (27 in four and six dimensions). We were able to do so 
for quite a few, but there are a still a number of codes that we were unable 
to do this for. In general we attempted this calculation using 500 decimal 
places of the maximum inner product. 

Moreover, for the minimal polynomials we calculated, we were able to 
verify all but two of them using the methods of Section HI These are 21 
points in five dimensions and 14 points in six dimensions, both of which have 
polynomials of high degree. The case of 27 points in four dimensions should 
also be noted because it has three rattlers, and since their positions are not 
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fixed relative to the other points we did not include them in our verification 
of the minimal polynomial. 

Lastly, we have a single new minimal polynomial in three dimensions, 
concerning the 22-point code. The minimal polynomial we found is 486u 18 + 
13113u 17 + 114996u 16 + 117476m 15 + 658256m 14 + 378752u 13 - 347056u 12 - 
121388m 11 + 81724m 10 - 70886u 9 - 55992u 8 + 12716m 7 + 6528m 6 - 2392m 5 - 
208m 4 + 284m 3 + 14u 2 + 5m + 4. But as in the two cases mentioned previously, 
we have been unable to verify this through finding an exact Gram matrix. 

Note that Tables [1] and [2] are the same codes as recorded in p~7] but 
with minimal polynomials associated with them. Also, some codes predate 
references to [5], but it is a convenient, systematic reference work. 



Table 1: Four-Dimensional Codes 



N 


U (20 decimal places) 


Minimal Polynomial 


Ref 


9 


0.16201519961163454918 


16m 3 - 16m 2 - Au + 1 


m 


10 


0.16666666666666666667 


6m- 1 





11 


0.23040556359455544174 


8m 3 - 12u A - 2u + 1 




12 


0.25000000000000000000 


4m - 1 





13 


0.30729565398102882233 


5632m 9 + 9472m* - 3072m 7 - 5888m 5 + 544m 5 + 
944m 4 + 152m 3 - 44m 2 - 14m - 1 




14 


0.31951859421260363550 


58m 7 + 174m 6 + 140m 5 - 16m 4 - 
54m 3 - 10m 2 + 4m + 1 




15 


0.35099217594534630330 


36m 4 - 18m 3 + 10m 2 - 1 




16 


0.38762817712253427776 


256m 1u + 1024m 9 + 256m 8 - 1152m 7 + 160m 6 + 
176m 5 - 144m 4 + 36m 3 + 37m 2 - 6m - 3 




17 


0.41226632322755925382 


Unknown 




18 


0.42281941407305934403 


14424m 16 + 42932m 15 + 18232m 14 - 62100m 13 - 
53831m 12 + 41528m 11 + 46442m 10 - 18248m 9 - 
20977m 8 + 6180m 7 + 5372m 6 - 1556m 5 - 
721m 4 + 240m 3 + 34m 2 - 16m + 1 




19 


0.43425854591066488219 


3m 2 + m - 1 




20 


0.43425854591066488219 


3m 2 + m - 1 


m 


21 


0.47138085850731791682 


16m 8 - 128m 7 - 64m 6 + 32m 5 + 
72m 4 + 32m 3 - 16m 2 - 8m + 1 




22 


0.49788413084355235629 


Unknown 




23 


0.50000000000000000000 


2m- 1 




24 


0.50000000000000000000 


2m- 1 


m 


25 


0.53731605665507787660 


Unknown 




26 


0.54078961769753707673 


3392m 6 + 2112m 5 - 496m 4 - 
656m 3 - 132m 2 + 6m - 1 




27 


0.55759135118017018253 


794m 5 + 393m 4 - 344m 3 - 82m 2 + 6m + 1 
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Table 2: Five-Dimensional Codes 



N 


u (20 decimal places) 


Minimal Polynomial 


Rcf 


11 


0.13285354259858991809 


45m 3 - 25u 2 - bu + 1 


i 


12 


0.15393160503302123095 


25m 4 + 30u 3 + 24m 2 +2u-l 




13 


0.18725985188285358702 


— — i n 

17u 3 - 5m 2 - 5m + 1 




14 


0.20000000000000000000 


5m- 1 




15 


0.20000000000000000000 


5u- 1 




16 


0.20000000000000000000 


5u- 1 


i 


17 


0.27047583526857362209 


~ ^ ~- 1 ~ T 

9m 4 - 16u 3 - 10u 2 + 1 




18 


0.27550174165981923839 


484m 5 - 488u 4 + 97u 3 + 17m 2 - u - 1 




19 


0.29182239902449014615 


57m 6 - 38m 5 - 109u 4 + 32m 3 + 23u 2 - lQu + 1 




20 


0.29938569289912478230 


5m 3 + 13u 2 -u-1 




21 


0.31491695717530346285 


869312m 14 + 8798656m 13 - 1062776m 12 - 
10586775W 11 - 968269m 10 + 3532907m 9 + 
188249m 8 - 659974m 7 - 11746m 6 + 72246m 5 - 
806m 4 - 5267m 3 - 97m 2 + 207u + 21 




22 


0.35499503416625620683 


Unknown 




23 


0.36977269694307633377 


Unknown 




24 


0.37423298246516725173 


1620m 1u + 5508u y - 5751m 8 - 2406u Y + 2055w b + 
276m 5 + 559m 4 - 210u 3 - 127u 2 + 48m - 4 
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Table 3: Six-Dimensional Codes 



N 


u (20 decimal places) 


Minimal Polynomial 


Rcf 


13 


0.11307975214744721385 


96m 3 - 36m 2 - 6u + 1 


m 


14 


0.13249092032347031437 


237299200m 13 + 463738880m 12 + 366062080m 11 + 
142462784m 10 + 23021376m 9 - 2398832m 8 - 
1849600m 7 - 343428m 6 - 1802m 5 + 11443m 4 + 
1390m 3 - 128m 2 - 20m + 1 




15 


0.14494897427831780982 


20m 2 + 4m - 1 




16 


0.17114764942939365334 


4000000m* + 15120000m 7 + 9896400M b - 
2451600m 5 - 4503600m 4 - 1511460m 3 - 
87480m 2 + 43740m + 6561 




17 


0.18327433702314481858 


400m 4 + 240m 3 + 16m 2 - 8m - 1 




18 


0.19781218860197545241 


784m 8 - 280m y - 119lM b + 120M b + 
218m 4 - 62m 3 - 31m 2 + 2m + 1 




19 


0.20022602589120548304 


Unknown 




20 


0.21428571428571428571 


14m -3 




21 


0.24285284801369170251 


Unknown 




22 


0.24886966078882474478 


Unknown 




23 


0.25000000000000000000 


4m- 1 




24 


0.25000000000000000000 


4m- 1 




25 


0.25000000000000000000 


4m- 1 




26 


0.25000000000000000000 


4m- 1 




27 


0.25000000000000000000 


4m- 1 
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Table 4: Seven-Dimensional Codes 



N 


u (20 decimal places) 


Minimal Polynomial 


Rcf 


15 


0.09870177627236447933 


175m 3 - 49m 2 -7u+l 





16 


0.11332087960014474125 


85823999m 12 + 153503984m 11 + 163426326m 1() + 
92426704m 9 + 29708081m 8 + 9292960m 7 + 
2252404m 6 + 410976m 5 + 29489m 4 - 
3216m 3 - 746m 2 - 48m - 1 




17 


0.12484758986639552862 


Unknown 




18 


0.12613198362288317392 


47m 2 + 2m - 1 




19 


0.15659738541709551030 


1280M e - 352m 4 - 48m 3 + 37m 2 + 3m - 1 




20 


0.16952084719853722593 


23m 2 + 2m - 1 


m 


21 


0.18152396080041583541 


Unknown 




22 


0.18274399763155681015 


19m 2 + 2m - 1 




23 


0.18274399763155681015 


19m 2 + 2m - 1 




24 


0.18274399763155681015 


19m 2 + 2m - 1 





Table 5: Eight-Dimensional Codes 



N 


m (20 decimal places) 


Minimal Polynomial 


Rcf 


17 


0.08773346332333854567 


288m 3 - 64m 2 - 8m + 1 


m 


18 


0.09946957270878709386 


Unknown 




19 


0.11140997502603998543 


Unknown 




20 


0.11949686668719356518 


244m 3 + 60m 2 - 2m - 1 




21 


0.13060193748187072126 


28m 2 + 4m - 1 




22 


0.13060193748187072126 


28m 2 + 4m - 1 




23 


0.15745541772761612286 


Unknown 




24 


0.15769214493936087799 


4m 4 - 4m 3 - 27m 2 - 2m + 1 
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